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Abstract. Cyg X-l exhibits irregular X-ray variability on all 
measured timescales. The usually applied shot noise models 
describe the typical short-term behavior of this source by su- 
perposition of randomly occuring shots with a distribution of 
shot durations. We have reanalyzed EXOSAT ME observations 
of Cyg X-l using the more general Linear State Space Models 
(LSSMs). These models, which explicitly take the observation- 
noise into account, model the intrinsic system variability with 
an autoregressive (AR) process. Our fits show that an AR 
process of first order can reproduce the system variability of 
Cyg X-l. A possible interpretation is again the superposition 
of individual shots, but with a single relaxation time r. This 
parameter was found to be 0.19 s. 

Key words: X-rays: stars - binaries: close - Accretion, accre- 
tion disks - Stars: Cyg X-l - Methods: statistical 



1. Introduction 

The X-ray source Cyg X-l was discovered by a rocket flight in 
1964 ( |Bowyer et al. 1965[ ). In March 1971 the sudden appear- 
ance of an accurately locateable radio source coinciding with a 
change in the luminosity of the X-ray source led to its identi- 
fication with the 09.7 lab supergiant HDE 226868 (Hjellming 
& Wade 1971; Bolton 1972). This star is known as part of a 
single lined spectroscopic binary system with an orbital period 
of 5.6 days at a distance of about 2.5 kpc. The X-ray emission 
of Cyg X-l is produced by the accretion of material from the 
supergiant primary onto the compact object. According to Her- 
rero et al. (1995), the most probable masses of the binary part- 
ners are 18 M Q for the primary and 10 M Q for the compact 
object, which is one of the most firmly established black hole 
candidates (BHCs). The accretion process is believed to be fu- 
eled by a focused stellar wind from HDE 226868 (Friend & 
Castor 1982). 

During most of the time Cyg X-l is emitting a non-thermal 
or hard state spectrum (for the soft state, cf. |Cui et al. 1997 ), 



which can be described by a power-law with a photon-index 
r « 1.5- 1.7, modified by an exponential cutoff with a folding 
energy of about 150 keV and reprocessing features. This spec- 
tral form can b e interpreted as b eing due to an accretion disk 
corona (ADC) ( Dove et al. 1997 and references therein). 

The short-term variability of the X-ray emission has been 
studied in order to gain better insight into the physical pro- 
cesses at work near the compact object. Cyg X-l was the first 
source for which X-ray variability on timescales < 1 s was de- 
tected (Odaet al. 1971). For along time this variability was sus- 



pected to be a special black hole signature, but today we know 
that X-ray binaries containing a ne utron star instead o f a black 
hole can display similar behavior ( Stella et al. 1985 ). Never- 
theless, the efforts in trying to identi fy a BHC by its irr egular 
short-term variability are continuing ( van der Klis 1995 ). 

Until now, no special process could be determined that de- 
scribes all the properties of the observed short-term variability. 
We define "variability on short timescales" as variability faster 
than a few 100 s, with emphasis on timescales < 1 s. These rapid 
fluctuations are usually described in terms of shot noise mod- 
els. It has become clear, however, that in the framework of ap- 
plying conventional shot noise models, complex shot profiles 
or distributions of shot durations an d amplitudes have to be as- 
sumed to explain the observations (Nolan et al. 1981 ; Lochner 
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et al. 1991; Negoro et al. 1994 ). Shot noise processes exhibit 
a fixed dynamical behavior in the sense that if the shots are 
given, no new dynamical information is produced during the 
run of the Poisson distributed point process. Only by the sum- 
mation of individual shots the temporal c orre lations in the time 
series of the process are generated (Sect. 2.3). 

We use the alternative method of applying Linear State 
Space Models (LSSMs) which are base d on stochast ic pro- 
cesses (i.e. autoregressive (AR) processes, Scargle 1981 ) to de- 
scribe the time series variability. In this case the dynamics of 
the process are produced using a different approach: each value 
of the time series refers to earlier values, with their temporal 
correlation being det ermined by the dynamical parameters of 
the system (Sect. 3.1). 

The time series of a standard shot noise (exponential decay 
with one relaxation timescale) and an AR process of first order 
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might look very similar, although the underlying processes dif- 
fer. The similarity in the time domain leads to shot noise and 
AR processes both exhibiting an exponentially decaying auto- 



Table 1. EXOSAT ME observations that have been analyzed 
with LSSMs. 



correlation function ( Konig & Timmer 1997). In contrast to the 



similarity of the time series, however, the theoretical frequency 
spectra are different (cf. Eq. ^|and|^). Whereas shot noise pro- 
files with varying relaxation timescales have to be added to 
approximate the observed periodogram ( |Lochneret al. 1991[ ), 
only one dynam ical AR parameter is needed to reproduce its 
properties (Sect. 3_3 and |4.l| ). 

We have studied the high time resolution EXOSAT ME 
lightcurves of Cyg X-l. The data are presented in Sect. 2.1. A 



commonly used method to analyze X-ray variability is to work 
in the frequency domain by fitting theoretical power spectra to 
the periodogram of a lightcurve: The periodogram of Cyg X-l 
and the power spectra of the usually applied shot noise mod- 
els are reviewed in Sect. 2.2 and 2.3, respectively. Our analysis 
of the lightcurves with LSSMs, however, is not based on fit- 
ting the periodogram but on an alternative procedure working 
in the time domain. Sect. |3| deals with the theory of the LSSM 
analysis. Our results are presented in Sect. |], and in Sect. ^ a 
possible explanation for the derived relaxation timescale is dis- 
cussed combining the temporal LSSM results and simulations 
of comptonized radiation using an accretion disk corona model 
for Cyg X-l. 



2. The short-term X-ray variability of Cyg X-l 

2.1. The data 

The EXOSAT raw data have been stored on Final Observation 
Tapes (FOTs) and are now available at the HEASARC archive. 
Table [l] lists the FOT observations of Cyg X-l that we have ana- 
lyzed with Linear State Space Models. We have chosen the ME 
datastreams provided by the primary timing telemetry modes 
HTR3 und HER6. These observations only contain events reg- 
istered in the Argon counters (1-20 keV). The lightcurves are 
given as countrates normalized to one half of the detector ar- 
ray (i.e. four Argon counters). Using the Interactive Analysis 
(IA) software, we extracted lightcurves corrected for dead-time 
effe cts and collimator efficiency (for an overview of the IA 
see 



Parmar et al. 1995 ). For the purpose of this paper an ex- 
plicit background subtraction is not necessary since the LSSM 
is implicitly modeling the measurement process (see Sect. 3.2, 
Eq.0). 

All EXOSAT observations of Cyg X-l have found the 
source in its usual hard state. An example for the character- 
istic hard state variability of Cyg X-l on short timescales can 
be seen in Fig. |l| 

2.2. The periodogram of Cyg X-l 

The periodogram P(fk) of a lightcurve y(tj) (where j = 
1, . . . , N) is providing the "strength" of harmonic variations 
with a certain frequency ff. in the lightcurve. P(fk) is de- 
fined as the squared modulus of the discrete Fourier transform 
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resolution [ms] 
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EE8087(05)* 


28.07.83 22:51 


07.813 
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EE1723(06)* 


25.09.83 05:15 


01.953 
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EE2461(06) 


21.05.84 04:59 


07.813 
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EE2956(06) 


07.07.84 10:57 


15.630 
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EE2956(07) 


07.07.84 13:23 


15.630 
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EE2929(10) 


09.07.84 14:59 


11.720 
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EX3620(08)* 


24.07.84 20:20 


07.813 
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EE3962(13) 


02.11.84 19:20 


09.766 
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EE3962(16) 


02.11.84 20:27 


11.720 



a Final observation tape number and number of observation therein 
(*: telemetry mode HTR3, else: HER6). 
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time [sec] 

Fig. 1. a X-ray lightcurve of Cyg X-l, showing a part of the 
fifth observation on FOT Nr. EE8087 (No. 1 of Table |j} energy 
range=l-20keV, bin time=62.4ms, length=156 s). b Blow-up 
of the marked segment of subfigure a with a length of 3.9 s and 
the original bin time of 7.8 ms, showing a typical shot structure. 
Both lightcurves are plotted without error bars for reasons of 
clarity. 

of (y(tj) — y) and is calculated for the Fourier frequencies 
f k = k/(N ■ At) with bin tim e At and k g {1, . . . , M}, where 
M is the integer part of N/2 ( |Scargle 1982b : 
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Fig. |2] shows the logarithmically plotted sample peri- 
odogram (solid line) of the observation No. 1 (Table [I]). It has 
been obtained by averaging over periodograms P(fk) from 48 
different lightcurve segments (details see caption of Fig. fy. 
Since P(fk) is a \% -distributed random variable, its standard 
deviation is equal to its mean (van der Klis 1989). Therefore, 
an individual periodogram exhibits large fluctuations (Fig. 0, 
dots). Calculating the sample periodogram significantly re- 
duces the scatter and allows the possibility to estimate the the- 
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Fig. 2. Sample periodogram (solid line) of Cyg X-l, averaged 
over 48 individual periodograms, which have been obtained by 
dividing the original lightcurve (No. 1, Table [j]) into 48 seg- 
ments of 156 s each and calculating P{jk) for each segment 
(Eq. [l]). The values of a typical individual periodogram, fluctu- 
ating wildly, are also displayed (dots). 

oretical spectrum of the process responsible for the observed 
variability. 

The periodogram of the short-term variability of Cyg X-l 
in its hard state is well known and exhibits the following distri- 
bution of timescales (Fig. |2|): 

- For frequencies above 10 Hz white noise dominates the pe- 
riodogram which means that the variability of the light- 
curves on all corresponding timescales is almost equally 
strong. These fluctuations (i.e. photon statistics, particle 
background) are introduced by the measurement procedure. 
It is common practice to subtract a constant (corrected for 
deadtime effects) whic h represents this noise com ponent 
from the periodogram ( Eelloni & Hasinger 1990b| ; Zhang 
et al. 1995). 

- Towards lower frequencies the power of the variability in- 
creases within the frequency range of roughly 0.04-10 Hz. 
This behavior is called red noise. It is often modeled by a 
power law, f~ a . Variability of the /~ a -type is known from 
many X-ray binaries (van der Klis 1995) and also from ac 
tive galactic nuclei ( |McHardy 1989|). In the case of Cyg X 



1, a is approximately 1 (Nolan et al. 
1991). 
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Fig. 3. a Shot noise power spectra for three different values of 
r, calculated according to Eq. ^ (A=l, h,Q=l). b LSSM power 
spectra for the same values of r as in (a), calculated according 
to Eq. ^| (of =1). All power spectra are given in arbitrary units 
and At has been normalized to 1 . In subfigure b, the variance 
of the observation-noise, a^, has been set to to allow direct 
comparison of a and b. 

stabilities in the mass transfer process. Sometimes the low 
frequency noise is associated w i th absorption dips i n the 
lightcurve ( |Kitamoto et al. 1984| ; |Angelini et al. 1994| ). 
- Several authors reported a transient quasi-periodic oscil- 
lation (QPO) feature in the power spectrum with a cen- 



tral frequency of about 0.04 Hz (Kouveliotou et al. 1992 



Angelini et al. 1994] ; [Vikhlinin et al. 1994| ; Borkus et al. 
1995), whereas other authors have not found any evidence 
for QPOs ( [Belloni & Hasinger 1990b| ; Miyamoto et al. 
1992). Thus the significance of this feature is not yet clear. 
In Fig. ||no QPO feature is present. 



2.3. The shot noise model 

Terrell (1972) proposed that the structures observed in the 
lightcurves of Cyg X-l might be due to a shot noise process 
s(t), i.e. the superposition of randomly occuring shots with the 
shot profile h(t'): 



1981 



Lochner et al. S W = X) h ^ ~ <SH ^ 



(2) 



- An important feature indicating the stationarity of the short- 
term variability process is the flattening of the power spec- 
trum for frequencies below «0.04Hz. This fiat top corre- 
sponds to the absence of additional long-term variations 
in the lightcurves. There are, however, some EXOSAT ME 
observations of Cyg X-l which show low frequency noise 
in the form of increasing power towards lower frequencies 



Here the tsH are the times at which the shots occur, with the 
time intervals between the isH following a Poisson distribution. 

The standard shot profile /i(i')sT is identical for all shots 
and consists of an instantaneous rise to height ho and an expo- 
nential decay with the decay time r: 



below 0.001 Hz (Angelini et al. 1994). Since this compo 



MOST = 



nent is not always present, it is most likely not produced 
in the same physical process as the short-term variability. 
Suggested causes for this low frequency variabilities are in- 



, t' < 

h ■ e~*> ,t' > 



(3) 



Shot noise models are often evaluated by comparing their 
theoretical power spectrum to the observed power spectrum 
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dBelloni & Hasinger 1990b| ; |Lochneret al. 199l[ ). The theoreti- 
cal power spectrum of the standard shot noise process is (Lehto 
1989): 



5*(/)sn oc 



(1/r) 2 + {2i: ff 



(4) 



where 1/A is the average time interval between the shots. 
S(/)sn provides the fiat top for / <C 1/(2itt) but has a fixed 
logarithmic slope of a —2 for higher frequencies (Fig. ||a). 
Since the shot noise process s(t) is continuous in i, whereas 
the observed lightcurve is the result of integrating the measured 
counts over a finite number of intervals At, S(ui)s~n must be 
corrected for binning. ^ochner et al. (1991)| have shown, that 
binning makes the shot noise power spectrum even steeper for 
high frequencies. 

Neither the observed slope nor the white noise floor of the 
Cyg X-l periodogram (Fig. ^) can be reproduced by the stan- 
dard shot noise process. In frequency domain fits the white 
noise level is usually treated as an additional constant parame- 
ter (cf. Sect. 22). To model the observed slope, different shot 
profiles with distributions of shot durations and shot ampli- 
tudes have been proposed. A number of those models, each 
having many degrees of freedom, are able to reproduce the ob- 
served periodogram and other second or der statistics (Belloni 
& Hasinger 1990b, Lochner et al. 1991 ). Higher order statis- 



tics like t he time skewness ( Priedhorsky et al. 1979 ) or phase 
portraits ( Lochner et al. 1991 ) were also studied but no special 
shot noise model could be singled out that allows for a homo- 
geneous description of different observations. In Sect. 3.3 we 



discuss the theoretical power spectrum of a first order LSSM 
and show that it can reproduce all the features of the hard state 
periodogram of Cyg X-l, requiring only one temporal parame- 
ter t (Fig. |b). 

3. The Linear State Space Model (LSSM) 

The LSSM analysis is a recently developed tool to model 
stochastic time series ( pantert 1993 ). It has e.g. been used to 
describe medical time series (hand-tremor data, Timmer 1995). 



Konig & Timmer (1997) were the first to apply the new method 



to astronomical data by successfully fitting lightcurves of active 
galactic nuclei with first order LSSMs (see also Konig et al. 
1997). The mathematical background of the LSSMs and the 
associated fitting-proced ure have been discussed in detail by 
Konig & Timmer (1997) and in the references therein. Here, we 
shall concentrate on first order LSSMs, for we found them to be 
appropriate to describe the Cyg X-l lightcurves (see Sect. Q). 

3.1. The autoregressive process 

One possibility to model irregular variability (i.e. temporal 
variations that cannot be forecasted exactly) is to assume that 



it is caused by a chaotic process ( Voges et al. 1987; Unno et al 



1990). The LSSMs, on the other hand, are based on the al- 
ternative assumption, that the variability of the observed sys- 
tem is produced by a stochastic process. The LSSMs use the 



rather general autoregressive (AR) formulation for the stochas- 
tic system-lightcurve x(tj). 



The AR processes were first introduced by Yule (1927) to 



model the variability of Wolf's s unspot numbe rs and have been 
well studied since then (e.g. by Scargle 1981 ). An AR process 
of order p (AR[p]) is defined by: 



x(t 3 ) 



i=l 



■ x(tj -i-At) + efo), e(t) oc A/"(0, a\) (5) 



The time series x(tj) is sampled at discrete times tj with 
time resolution At. The purely stochastic component e(t), is a 
Gaussian random variable with mean and variance of . Since 
Eq. (Q) is of the same structure as a regression equation for the 
variables x(tj) and (x(tj — i ■ At)) the name autoregressive 
process has been assigned to it. 

In an AR lightcurve, for each tj the value of the random 
variable x(tj ) is correlated with the values of the process at ear- 
lier times. This correlation is decreasing with increasing time 
differences between two lightcurve values (expressed by an ex- 
ponentially decaying autocorrelation function). The stochastic 
component e(t) is the reason that the process does not sim- 
ply come to rest. e(t) is an intrinsic property of the system- 
variability: the system-noise. 

The temporal correlations in the AR lightcurve are char- 
acterized by the dynamical parameters (a^). The p dynamical 
parameters are related to p temporal parameters, describing the 
temporal structures in the lightcurves: depending on the pro- 
cess, the (cii) represent stochastic relaxators with relaxation 
times r, or damped stochast ic oscillators with relaxation times 
r and periods P, or both ( Honerkamp 1994 ). In case of an 
first order AR process there is only one dynamical parameter 
ai = a. For stationary processes \a\ must be <1 and only pos- 
itive values of a lead to plausible physical models. The corre- 
sponding temporal parameter, t = — 1/ In \a\, is the relaxation 
time of a stochastic relaxator, one representation of which is an 
exponentially decaying shot as described by Eq. (||). 



3.2. The observation-noise 

For the determination of the temporal parameters of the system- 
process the noise which is caused by the measurement (i.e. pho- 
ton statistic, particle background) has to be considered because 
it disturbs the temporal structures of the system-lightcurve. 
If the observation of a system-lightcurve is modeled with a 
plain AR process, the so called observation-noise will lead 
to an underestimation of the true temporal system-parameters 



(Robinson & Nather 1979; Konig & Timmer 1997). In order to 



solve this general problem, a LSSM consists of two equations: 
the system-equation for the intrinsic system-lightcurve and the 
observation-equation for the measured lightcurve. The latter 
describes the observed lightcurve y(tj) by explicitly consider- 
ing the influence of the observation-noise r)(tj) on the system- 
lightcurve. 
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The important model for the analysis of the short-term vari- 
ability of Cyg X-l is the LSSM of first order. 

x(tj) = a ■ x 

y{tj) = c ■ x(tj) + vitj), nitj) oc Af(0, o*) (7) 



system-lightcurve x[tj) by using the observed lightcurve y(tj) 
and a given set of LSSM parameters. In the following maxi- 
mization step the likelihood function is maximized, taking the 
estimated x(tj) into account and providing a new set of pa- 
rameters. The results of this iterative procedure are best esti- 
mates for the LSSM[p] parameters and for the intrinsic system- 



Like the system-noise e{t,\ the observation-noise jjfo) is lightcurve ( a detailed example is discussed in Sect. O). Note, 



modeled by a Gaussian noise component; c is a constant nor 
malization factor. The equations for higher order LSSMs are 
given by Konig & Timmer (1997) (the temporal parameters of 
an LSSM of order p (LSSM[p]) correspond to the p dynamical 
parameters of Eq. |[). 

3.3. The power spectrum 

The power spectrum S'(/)lssm[i] of a first order Linear Stat e 
Space Model has the following form ( Konig & Timmer 1997| ): 



>5'(/)lssm[i] k 



— 2 acos(27r/) 



(8) 



The examples for S'(/)lssm[i] displayed in Fig. ||b show the 
flat top as well as red noise. The observational white noise 
floor is also provided by Eq. || but has been omitted in Fig. |]b 
(ct^ =0) to allow the direct comparison with the standard 
shot noise spectra (Fig. |3|a): In contrary to the fixed red noise 
slope a =2 of the standard shot noise spectra, the slope of 
>5(/)lssm[i] can be modeled by adjusting a, i.e. r. All the fea- 
tures of the hard state power spectrum of Cyg X-l can thus 
be reproduced by the power spectrum of the LSSM[1] which 
is defined by only one temporal parameter r = — 1/ log |a|. 
Furthermore, Eq. [8] does not need to be corrected for binning 
effects, since the definition of the LSSM already is discrete in 
t (see Eq. ||and0). 

Principally, LSSMs can be evaluated by fitting their power 
spectra to the measured periodogram. However, fits in the fre- 
quency domain do not allow the explicit modeling of the obser- 
vational noise and may be influenced by several other sources 
of uncertainties, especially by spectral leakage, which contam- 
inates the power spectrum of each finite time series (Deeter & 
Boynton 1982; Deeter 1984). To avoid these uncertainties our 
LSSM analysis was performed in the time domain by direct 
comparison of the model-lightcurves and the measured light- 
curves. This method allows fits with a statistical significance 
higher than that of the frequency domain fits (Konig & Timmer 
1997). 



3.4. The fitting procedure 

The parameters of a LSSM[p] fit are estimated with the help 
of the maximum likelihood procedure: the set of parameters is 
derived, for which the probability of observing the measured 
lightcurve y(t) is at maximum. We have obtained those pa- 
rameters using the expectation-maximization (EM) algorithm 



(Honerkamp 1994). The method is iterative, starting from a set 



that the latter is not produced using a smoothing filter: all vari- 
ability timescales of the observed lightcurve (see e.g. Fig. ^) 
are still present in the derived system-lightcurve (see e.g. 
Fig- 0b). 

The length and the time resolution of lightcurves, which are 
analyzed by LSSMs, have to be chosen considering the follow- 
ing requirements: 

- The duration of the lightcurve has to be several times the 
relaxation timescale, which is expected to be a few tenths 
of a second (shot noise models), or the relaxator cannot be 
found. 

- The bin time has to be well below the relaxation timescale 
to allow us to distinguish the relaxator from statistical fluc- 
tuations. 



Simulations by Konig (1997) have shown that the duration 



should be at least 5 times larger and that the bin time should be 
at least 10 times smaller than r. In order to fulfill these criteria 
without an unreasonable increase in LSSM computing time, we 
have fitted lightcurves with a bin time of 16 ms and length of 
16 s (e.g. Fig. |]a). Although this time interval is shorter than the 
intervals typically studied in the context of shot noise analysis, 
it is long enough for the determination of the relaxator. 

4. Results 

4.1. LSSM fits of various order to one exemplary lightcurve 

Table § shows the results of LSSM fits of order 0-5 to the light- 
curve segment displayed in Fig. [|a: The relaxation timescale 
t of about 0.2 s, which is derived by the LSSM[1] fit, is also 
found by the fits of higher order (the LSSM[0] represents un- 
corrected white noise and therefore does not yield a temporal 
parameter). All other temporal correlations are decaying within 
a few time bins. Since they are due to statistical fluctuations 



they are negligible (Konig 1997). This means, that the dynam 



of initial parameter values. In the expectation step, the Kalman 
filter is applied, which allows to estimate the unobservable 



ics of the lightcurve can most likely be described by one re- 
laxation timescale. If this is the case, the LSSM[1] fit should 
deliver an adequate description of the observed lightcurve and 
the goodness of fit should not increase significantly for higher 
order LSSMs. 

A statistical test on the residuals, i.e. the difference between 
the derived system-lightcurve and the observed lightcurve, was 
performed to determine the quality of the LSSM fits: If a given 
LSSM can describe the observation, the residuals should be an 
uncorrelated white noise realization (Eq. [7]). In this case all the 
temporal correlations of the observed lightcurve are modeled 
by the system-lightcurve. We use the Kolmogorov-Smirnov 
(KS) test to quantify the goodness of the LSSM fits. It com- 
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Fig. 4. a One of the 900 EXOSAT ME Cyg X-l lightcurves 
with a length of 16 s and a bin time of 16 ms, which were mod- 
eled by LSSMs. Mean countrate: 633.1 counts/sec/half, stan- 
dard deviation: 234.5 counts/sec/half, b Best estimate for the 
autoregressive system-lightcurve, which has been obtained by 
applying a LSSM of order p =1 to the lightcurve displayed 
above (see Table Both lightcurves are shown without error 
bars for reasons of clarity. 

Table 2. Results of LSSM fits of various order to the lightcurve 
in Fig. 



Order p of 
the LSSM 


Relaxator r 
M 


Period P 

M 


KS probability 

% 









0.0 


1 


0.185 




93.5 


2 


0.209 
0.018 




87.0 


3 


0.177 
0.022 


0.092 


97.6 


4 


0.204 
0.046 
0.011 


0.128 


88.5 


5 


0.175 
0.046 
0.015 


0.128 
0.096 


97.4 



pares the power spectrum of the residuals to a flat white noise 
power spectrum. 

In Table ^| we list the KS probabilities for the residual light- 
curves to actually be white noise realizations. For the LSSM[1] 
fit a very high probablity of 93.5% is reached, which means that 
the observed lightcurve can be well modelled by an LSSM[1] 
process. This is consistent with the visual impression from 
Fig. the LSSMfl] estimate (b) tracks the observed lightcurve 
(a) very closely. Higher order LSSM fits do not improve the 
the fit significantly. The LSSM[0] (white noise) is of course 
rejected at any level of confidence. 

The distribution of the LSSMfl] model residuals is shown 
in Fig. |5|a. Its mean and variance are 0.7 counts/sec/half and 
(179.9counts/sec/half) 2 , respectively. Thus the variance of the 
observation-noise is about 82% of the upper limit as given by 



40 



! (a) 

i aA. 








- (b) 


' * * J 



-1000 -500 500 1000 

residuals [counts/sec/half] 



Fig. 5. a Distribution and b normal quantile plot of the residuals 
of the LSSM[1] fit to the lightcurve given in Fig. ^a. The lines 
in subfigure a indicate the mean and standard deviation of the 
observation-noise. The straight line in subfigure b corresponds 
to a Gaussian distribution having the same mean and variance 
as the distribution of the residual lightcurve. 

the mean countrate of the original lightcurve and assuming a 
pure Poisson process. The lower half of Fig. || displays the nor- 
mal quantile plot of the fit residuals: it arranges the residuals 
in increasing order and for each data value indicates the posi- 
tion that corresponds to the same probability in a normal dis- 
tribution. If the data are normally distributed, all points should 
lie on a straight line. Fig. |5|b shows that the distribution of the 
residuals can be approximated by a Gaussian distribution, but is 
compressed for low and stretched for high countrates. This in- 
dicates that the fit residuals are Poisson-distributed as expected 
for the observation-noise (mean countrate of the original light- 
curve: 10.1 counts per 16 ms timebin). The LSSM[1] model 
parameters reflect the Poisson distribution of the observation- 
noise by underestimating the Gaussian variance cr 2 given as 
(143.8 counts/sec/half) 2 in this example. 

4.2. LSSM fits of first order to a sample of lightcurves 

In order to obtain a reliable value for the temporal parame- 
ter t, we have fitted LSSM[l]s to a sample of 900 lightcurve 
segments covering a total length of four hours. We have used 
1600 s of each of the nine observations listed in Table [j]. The 
1600 s consist of 100 equally long, uninterrupted segments with 
a temporal resolution of 16 ms. 

For each of the nine observations 100 values for the relax- 
ator t were obtained. In Table ||, the mean f, the standard de- 
viation s T , and the most probable value r max of these distribu- 
tions are listed. Averaging f , s T , and r max over the nine obser- 
vations gives (0.28±0.04)s, (0.13±0.03)s, and (0.19±0.03)s, 
respectively. All distributions exhibit the same asymmetric 
form, in the sense that r max is always lower than f but still 
higher than the lrr deviation. This kind of distribution is consis- 
tent with the results that have been obtained from simulations 
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Table 3. Results of the LSSM(AR[1]) fits to EXOSAT ME ob- 
servations of Cyg X- 1 . 



No. a energy channels 



count- 
rate 13 



T 
M 



St Tmax 

[s] [s] 



1 


1-128 


565±57 


0.26 


0.10 


0.22 


2 


1-128 


567±37 


0.36 


0.17 


0.22 


3 


9-21, 22-51 


541±33 


0.30 


0.16 


0.22 


4 


5-19, 27-65 


388±19 


0.26 


0.12 


0.17 


5 


5-16, 27-65 


457±29 


0.31 


0.16 


0.19 


6 


7-21, 26-81 


499±35 


0.25 


0.11 


0.17 


7 


1-128 


389±24 


0.25 


0.09 


0.22 


8 


5-16, 17^9 


359±14 


0.27 


0.14 


0.17 


9 


5-16, 17^9 


321±14 


0.24 


0.11 


0.17 



erage over the mean countrates of the 100 lightcurve segments and its 
standard deviation ([counts/sec/half]). 



150 - 




50 - 



0.4 0.6 0.8 

i [s], bin width = 0.048 s 



Fig. 6. Distribution of values of the relaxation timescale r, ob- 
tained from LSSM[1] fits to a sample of 900 Cyg X-l light- 
curves. The mean value of the distribution is f =0.28 s (dash- 
dotted line) with the standard deviation s T =0.14s (dashed 
lines). The most probable value is r max =0.22 s. 

of AR lightcurves based on one temporal parameter r (Hamil- 
ton & Wu 1987). 

The lightcurves of three of the observations (No. 1, 2, and 
7 of Table |3]) take the whole Argon energy range (l-20keV) 
into account, i.e. their absolute countrates may be compared. 
Even though the analyzed segments of one of those three ob- 
servations (No. 7) have a considerably smaller mean countrate 
than those of the other two, the corresponding distribution of 
t is not significantly different. This means that no correlation 
between t and the luminosity can be seen, which is confirmed 
within each of the nine observations: the mean countrates and 
the relaxation times of the 100 segments are not correlated. 

Since the observations provide consistent results, all 900 
values for r can be combined, in order to arrive at a distribu- 
tion with better statistics (Fig. ^). Note, that this non-Gaussian 
distribution, resulting from the maxim um likelihood estimatio n 
of r, cannot be described analytically ( Hamilton & Wu 1987 ). 
In addition, simulations have shown that the distribution of r is 



contaminated by bad fit results at its "high r end" due to numer- 
ical effects ( Hamilton & Wu 1987[ |Konig 1997[ ). Therefore, f 
can be expected to overestimate the "true" relaxation timescale 
of Cyg X-l, whereas the most probable value T max should be a 
better estimator. 



5. Summary and Discussion 

5.7. The LSSM[1 ] results — a generalization of the shot noise 
approach 

We have shown that the short-term variability of the EXOSAT 
ME lightcurves of Cyg X-l can be well described by a Lin- 
ear State Space Model of first order, i.e. the dynamics of the 
system can be modeled as an autoregressive process with one 
temporal parameter - the relaxation timescale r. We find r to 
be (0.19±0.03)s for the ME energy-range. The relaxator is not 
found to be correlated either with time nor with the luminosity 
of the source, indicating that the physical process producing the 
emitted radiation is very stable and suggesting that the short- 
term variability is independent of the mass accretion rate of the 
system. The LSSM analysis in the time domain allows us to 
estimate the fraction of the countrate that is due to observation- 
noise. When analyzing the periodogram in the frequency do- 
main this white noise component can only be represented by 
a constant whose extraction is not trivial (Belloni & Hasinger 
1990a; Zhang et al. 1995). Furthermore, the problem of spec- 
tral leakage that has to be delt with in the frequency domain, is 
circumvented by fitting in the time domain. For these reasons 
our LSSM analysis delivers results with higher statistical sig- 
nificance than corresponding frequency domain fits (Konig & 
Timmer 1997). 

These LSSM[1] results allow a much simpler description 
of hard state short-term variability than multi-timescale shot 
noise models. Although a detailed quantitative comparison is 
beyond the scope of this paper, the reproduction of the pe- 
riodogram by the LSSM[1] (compared to its approximation 
by adding shot-profiles with different timescales, corrected for 
observation-noise and binning) as well as the greater sensitivity 
of the LSSM fitting procedure (compared to frequency domain 
fits, which are usually used to evaluate shot noise models) sug- 
gest that LSSM[l]s are better suited to describe the nature of 
the observed variability. 

The LSSMs can model the different realizations of a 
stochastic relaxator t, whereas shot noise models are generally 
restricted by the definition of special shot forms. Shot noise 
lightcurves therefore might be regarded as a subclass of LSSM 
lightcurves in the sense that a superposition of exponentially 
decaying shots can be interpreted as one possible realization 
of an intrinsic AR[1] process. The inspection of the measured 
lightcurves of Cyg X-l implies that the source of the derived 
AR[1] process is indeed the stochastic superposition of individ- 
ual shot events, corresponding to the basic idea of shot noise. 
We note that on larger timescales this kind of variability is also 
pres ent in the X-ray emission of active gala ctic nuclei (AGN) 
(e.g. |McHardy 1989| ; [Mushotzky et al. 1993] ; Konig & Timmer 
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1997; Konig et al. 1997). The physical mechanism responsible 
for such a temporal behavior, however, is not yet understood. 

5.2. The LSSM results in the light of time-dependent 
Comptonization models 

Recently, the discussion concerning accretion physics has be- 
gun to concentrate on the consideration of timing and spectral 
properti es of the X-ray emission as two aspect s of the same 
model (Kazanas et al. 1997; Konig et al. 1997; Wilms et al. 



1997; and references therein). The spectrum of both, AGN and 
the hard state of galactic black hole candidates, is usually ex- 
plained by inverse Comptonization, where soft X-ray photons, 
provided by a cold accretion disk, are upscattered by inverse 
Compton collisions in a hot plasma to produce the observed 
high energy power-law. 

In this context it can be assumed that the observed X-ray 
temporal behavior is the result of the processing of short shots 



of seed photons within the accretion disk corona (Payne 1980). 
The initial photon pulse is hardened and temporally broad- 
ened while diffusing through the corona. As harder photons 
in the emerging spectrum have, on average, undergone more 
scattering events, they have stayed in the corona for a longer 
time than softer photons and the eme rging pulse is com par- 
atively broader (observed in AGN by Konig et al. 1997 ). In 
addition, variability structures in high energy lightcurves ex- 
hibit a characteristic time-lag with respect to those in low en- 
ergy lightcurves (Fig. [7]). Frequency-dependent time-lags have 
e.g. been observed in Ginga and RXTE observations of Cyg X- 
1 ( [Miyamoto et al. 1988] ; [Miyamoto et al. 1992| ; Wilms et al. 
1997). The question to what extent these lags are produced in 
a hot corona or wether they are intrinsic properties of the cold 



disk is still unresolved ( 


Miyamoto et al. 1988 


Nowak & Vaughan 1996; 


Mowak et al. 1998). 



Miller 1995 



Nowak et al. 

(1998), however, show that the frequency-dependent time-lags 
of Cyg X-l can be qualitatively reproduced using a simple 
propagation model (but they also point out that the longest ob- 
served time-lags («0.1 sec) can probably not be explained by 
Comptonization). Another challenge for combined spectral and 
temporal theories is given by modelling the coherence of vari- 
ability structures in different energy bands which is observed to 
be near unity ov er a large frequency range in Cyg X- 1 (Vaughan 
& Nowak 1997; |Wilms et al. 1997| ). 

Additional support for the time-dependent Comptonization 
model comes from the joint spectral and temporal analysis of 
AGN, which have similar properties to those of Cyg X-l. For a 
sample of AGN observations we were able to show that a lin- 
ear correlation between the photon index T and the LSSM[1] 
timescale r exists, in the sense that harder spectra have a longer 
variability timescale (Konig et al. 1997). By comparing the 



observed T-t relation with Monte Carlo simulations of time- 
dependent Comptonization models, it is possible to scale the 



model geometry ( Konig 1997 ). For the hard-state galactic black 
hole candidates a similar study cannot be performed since not 
enough objects are known. What can be done, however, is to 
check whether our best-fit LSSM models for Cyg X-l could, at 



300 — 
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200 - 
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I 100 - 
50 - 



10 15 
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Fig. 7. Shot profiles for the energy-ranges 4-5 keV (dash- 
dotted line) and 12-23 keV (dashed line). The profiles have 
been computed for an accretion disk corona (ADC) model with 
an accretion geometry as described by pove et al. (1997) . This 
geometry is similar to an advection dominated flow (ADAF, 
Abramowicz et al. 1995| ; Marayan 1996). The ADC model pa- 
rameters were chosen to be appropriate for Cyg X-l (disk tem- 
perature: /cTd=200eV, coronal temperature: fcT c =66keV, opti- 
cal depth of the corona: r c =2.1; Dove et al. 1997). 



least roughly, be explained with a time-dependent Comptoniza- 
tion model. 

We have computed shot profiles for several possible ac- 
cretion geometries using a linear Monte Carlo code. The de- 
tailed results of these simulations will be discussed elsewhere. 
Since shot noise can be regar ded as a special representation 
of an AR[1] process (see Sect. 5.1), we can use the computed 
shots to generate lightcurves based on these profiles, includ- 
ing an additional white noise component. In analogy to our 
treatment of the distribution of r seen in the observations, we 
derive a most probable value for r from the simulated light- 
curve samples, Tth. We find Tth = (3 . . . 6)R c /c, with R c be- 
ing the radius of the corona. The profile of a Compton-shot 
only depends on the relative size of the system, parameter- 
ized by the light travel time through one radius of the spher- 
ical corona, R c /c. By identifying Tth with the measured value 
of r =0.19 s, it is therefore possible to express the measured 
r in terms of the coronal radius. Our simulations give a first 
estimate of 320-640 Schwarzschild radii for the coronal radius 
(assuming a mass of IQMq for the black hole). Other authors 
estimate the size of the Comptonization cloud in the hard state 
to be w 100 Schwarzschild radii (foin et al. 1997[ ADAFs) or 
«23 Schwarzschild radii (Nowak et al. 1998, from minimal 



time-lags). These values are only in rough agreement and a bet- 
ter understanding of the relation between spectral and temporal 
properties is needed to arrive at a consistent picture. We plan 
a more detailed study of the Compton-shot profiles as well as 
LSSM analyses of RXTE data to further constrain the accretion 
geometry of Cyg X-l. 
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